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ABSTRACT 

Optimal-control and finite-thrust aspects of the approach maneuver - planetocentric 
orbit - departure maneuver sequence at the destination planet of a round trip are studied 
Solutions are developed for single maneuvers, and conditions for joining these into an 
optimal total trajectory are derived. Numerical results are presented for an elliptic 
parking orbit and typical asymptotic velocities and directions. It is shown that "gravity 
losses" and "steering losses" are not directly additive, low vehicle thrust to weight ra- 
tios are desirable, and trajectories with single-burn escape and capture maneuvers are 
usually advantageous. 
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OPTIMAL FINITE-THRUST TRANSFER BETWEEN PLANET-APPROACH 
AND DEPARTURE ASYMPTOTES WITH SPECIFIED 


INTERMEDIATE ORBIT 
by Edward A. Willis, Jr. 

Lewis Research Center 

SUMMARY 

This report deals with finite-thrust and optimal -control aspects of the approach 
maneuver - planetocentric orbit - departure maneuver trajectory sequence at the desti- 
nation of an interplanetary round trip. Individual escape or capture maneuvers (with 
constant, continuous acceleration) leading from an arbitrary planetocentric orbit to a 
specified hyperbolic asymptote are first studied using the Pontryagin maximum principle. 
An analytical "matching condition" is then derived which permits an optimal approach - 
orbit - departure trajectory sequence to be constructed from single maneuver data. Rep- 
resentative results obtained in this manner are compared with alternative (multiburn) so- 
lutions to this same problem. It is shown that the present "single -burn" class of solu- 
tions is usually superior unless the asymptotic direction constraint is very severe. 

It is also found that "gravity losses" due to finite thrust and "steering losses" due to 
constraints on the asymptotic direction do not add together directly. In the medium to 
high acceleration regime, gravity losses occur primarily in connection with near- 
optimal-angle maneuvers (i.e. , with no asymptotic constraints). When the asymptotic 
direction is strongly constrained, there are large steering losses but little or no addi- 
tional penalty for finite thrust. It is concluded that the steering losses can be offset to 
a significant extent by using a relatively small, lightweight engine which need not develop 
vehicle accelerations larger than a few percent of local gravity. This is illustrated in 
terms of a specific mission (stopover round trip to Jupiter using nuclear rockets) . For 
the case considered, it is found that the most desirable vehicle thrust to weight ratio at 
Jupiter is in the range of 0.01 to 0.03 local gravity (0.02 to 0.06 Earth gravity). 


INTRODUCTION 


Many significant problems remain before an ambitious interplanetary round trip like 
the one in figure 1(a) can be accomplished. Trajectory problems in particular require 
early attention because their solutions comprise an essential "input" for hardware - 
oriented studies such as the selection of a desirable engine type and size for future de- 
velopment. One of the most interesting parts of the total trajectory is the multipoint 
boundary value problem arising at the destination planet of a round trip. There the space 
vehicle must transfer between prescribed hyperbolic asymptotes by way of an interme- 
diate parking orbit, as suggested by figure 1(b). The two asymptotes are defined at the 
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(a) Typical heliocentric trajectory (Jupiter mission). 


--Sphere of 
/ Influence 



(b) Pianetocentric 
trajectory at 
destination. 


Figure 1. - Elements of interplanetary trajectories. 


sphere of influence by the hyperbolic velocity vectors V OT ^ and ^ as shown in fig- 
ure 2(a). The intermediate orbit itself represents a third boundary condition which with 
one of the two asymptotes leads to the single -maneuver geometry illustrated in figure 2(b) . 
The problem, then, is to accomplish the indicated approach maneuver - intermediate 
orbit - departure maneuver sequence with minimum propulsive effort or AV. 

Solutions are easily obtained if the intermediate orbit is circular; the geometric 
boundary value aspect is then trivial, and the question of optimally escaping from a cir- 
cular orbit has already been discussed at length (e. g. , in refs. 1 to 4). But when an 
elliptic or hyperbolic orbit is used, the solutions are more interesting because several 
extra degrees of freedom are then relevant. Prior studies of the elliptic -parking- orbit 
case (refs. 5 to 7) consider only the impulsive thrust limit. These results (which are 
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applicable to stopover missions) indicate large AV reductions in particular examples 
for using an elliptic rather than a circular parking orbit. Even greater AV savings are 
obtained for nonstop round trips for which the intermediate orbit is hyperbolic. In this 
case also, prior studies such as references 8 and 9 are limited to the impulsive thrust 
approximation. These results are of undoubted value because in many cases they provide 
limits of performance and a qualitatively correct description of the problem. On the 
other hand, it is clear that more realistic methods are needed to support detailed mis- 
sion studies which might involve, for example, the determination of optimal engine sizes 
or vehicle thrust to weight ratios. 

Accordingly, the present analysis deals with finite-thrust and optimal-control aspects 
of the problem illustrated in figure 2. It is intended to be a direct extension of refer- 
ence 4, which treated only the optimal-angle case for individual (escape or capture) 
maneuvers. The question of minimizing AV is approached by separating the total tra- 
jectory into a pair of escape and capture maneuvers (as shown in fig. 2(b)). These are 
studied individually using Pontryagin's maximum principle (ref. 10) and then assembled 
by means of a matching condition to produce an optimal overall solution. 

For single-burn escape and capture maneuvers, the propulsive effort (AV) and the 
associated state and control variable histories are studied as functions of the magnitude 
and direction of V^. Limits of performance are defined in terms of closed-form varia- 
tional solutions valid in the limits of impulsive or very low thrust. (The impulsive solu- 
tion is also used as the basis of a stable and very efficient feedback algorithm; this fur- 
nishes a convenient means of approximating the variational results presented herein.) 

A simple analytical matching condition is derived from the principle of "equal incre- 
mental slopes". This permits the angular constraint to be distributed optimally between 
the escape and capture arcs (without the need for trial-and-error searching) and thereby 
yields an optimal overall trajectory. 

Illustrative numerical results are presented, and the global optimality of the present 
(single-burn) class of solutions is discussed by comparing them in an approximate man- 
ner with multiburn maneuvers which also yield low AV. 

Finally, an example is given (stopover round trip to Jupiter) to illustrate how the 
present results may affect important questions such as the selection of an advantageous 
engine type and size for future development. 


SYMBOLS 


A acceleration 

a dimensionless acceleration, see eq. (4) 
b auxiliary function used in appendixes A and C 
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Av 


X 


defined in eq. (20) 

terminal error function, see eq. (46) 
components of terminal error function, see eq. (46) 
eccentricity of conic section 
defined in eq. (11) 

gravity and steering loss AV correction factor, AV/AV m ^ n 
acceleration due to gravity 

dimensionless acceleration due to gravity, see eq. (4) 

angular momentum 

Hamiltonian function 

angular momentum/R^ V„ 

Pci c j pa 

specific impulse, sec 
thrust switching function 
mass 

mass fraction 

magnitude of primer vector, see eq. (16) 
semilatus rectum/Rp a 
polar radius 

dimensionless radius, see eq. (4) 

orbital period 

time 

thrust angle of attack 

thrust magnitude control function 

velocity 

propulsive velocity increment 

impulsive AV for unconstrained asymptotic direction, see eq. (32) 
dimensionless velocity, see eq. (4) 

propulsive velocity increment in dimensionless units, AV/V 

c, pa 

correlation angle, defined by eq. (54) 



Y auxiliary function defined in eq. (CIO) 

a trajectory path angle 

rj efficiency 

0 angle between hyperbolic asymptotes 

9 trajectory central angle measured from major axis of intermediate orbit 

ju planet gravitational constant 

v parking orbit true anomaly at thrust initiation (escape maneuver) or at thrust 

termination (capture) 

a mismatch angle, seep. 35 

r dimensionless time, measured from zero at the inner end of a trajectory posi- 

tively toward the sphere of influence 

cp true anomaly on escape or capture hyperbola 

i/a adjoint variables 

a) planet angular velocity, see fig. 2 

( )' denotes differentiation with respect to t 

( ) denotes mean value 

( )* denotes temporary normalization 

Subscripts : 

aa apoapse 

act actual 

c circular 

cst coasting 

des desired 

e engine 

g to be gained 

h junction between powered arc and approach or departure hyperbola 

i,j,k general numerical indices 

imp impulsive 

m mean 

max maximum 
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min 


minimum 


0 junction between powered arc and.planetocentric orbit 

opt optimal 

p propellant 

pa periapse 

pi planet 

po parking or intermediate orbit 

pr propulsive 

ref reference 

sc spacecraft 

°° sphere of influence 

% Jupiter 

® Earth 

1 capture maneuver at destination planet 

2 escape maneuver at destination planet 

ANALYSIS 

The stated problem involves trajectories which are consistent with the boundary con- 
ditions illustrated in figure 2 and for which the sum of the capture and escape propulsive 
velocity increments 


AV X + AV 2 

is a minimum. For example, the escape maneuver contributes an amount 


av 2 = 



Thrust ^ 
Mass 


( 1 ) 


to this criterion. This is directly related to space vehicle mass by the classical rocket 
equation which may be written as 
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( 2 ) 


M l,h = M 2,h ex P 


r Av i - av 2^ 


IG 


© 


when the propulsion system operates with constant jet velocity (i. e. , IG 0 ) and no mass is 
jettisoned between the escape and capture maneuvers. The approach used in minimizing 
equation (1) is, first, to study individual maneuvers which yield minimum AV for arbi- 
trary magnitude and direction of V w , and second, to develop matching conditions which, 
when satisfied jointly by the escape and capture maneuvers, will produce an optimal 
overall trajectory. 


Assumptions and Basic Equations 

The planetocentric orbit escape or capture maneuver as illustrated in figure 2(b) is 
studied under the following conventions : 

(1) In keeping with the usual "successive two-body" trajectory model, an interplanetary 
trip (fig. 1) is viewed as consisting of alternate heliocentric and planetocentric arcs 
matched at the sphere of influence of each planet as indicated in figure 2(a) . The radius 
of the sphere of influence is taken to be negligible in comparison with interplanetary dis- 
tances yet very much larger than a characteristic dimension of the plenetocentric orbit 
(for instance the apoapse of an ellipse) . This further implies that only the magnitude and 
direction of are specified. Its lateral displacement relative to planet center is taken 
as an open parameter for optimization. 

(2) The only external forces along the trajectory are the engine thrust and inverse- 
square gravitational attraction. 

(3) The entire planetocentric trajectory (fig. 1(b)) lies in a single plane passing through 
the center of force and containing the two asymptote vectors V ’ 1 and V w 

(4) The V 1 and V 2 vectors are regarded as being determined in magnitude and 
direction by prior heliocentric calculations such as those in reference 7. Thus, possible 
interactions between heliocentric and planetocentric trajectory optimizations are disre- 
garded here. 

(5) The previous four points further imply that the maneuver time (from initiation of 
the first burn until escape energy has been reached) must not exceed a small fraction of 
the time spent within the sphere of influence. If this were not so, then either the asymp- 
totes derived from interplanetary calculations would be invalidated, or a significant part 
of the planetary stay time would be spent in accomplishing propulsive maneuvers instead 
of performing mission objectives. 
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(6) Primary attention is directed toward the class of single-burn maneuvers; these are 
deemed preferable because of their simplicity and their short maneuver times. More- 
over, it will be demonstrated that they are usually superior to the multiburn solutions 
when both types are required to satisfy item (5). These single-burn trajectories consist 
of a continuous -thrust maneuver from the orbit departure point to burnout; a coasting arc 
then proceeds to the sphere of influence where it must yield the desired magnitude and di- 
rection of V . 

oo 

(7) Vehicle dynamics are neglected ; it is assumed that commanded control settings are 
attained without significant lag. 

(8) The continuously variable direction and bounded magnitude of the vehicle accelera- 
tion are taken as the primary control functions. Also to be optimized are the power-on 
and power-off points, the lateral displacements of the two V w vectors, and the orienta- 
tion of the intermediate orbit relative to these vectors. 

With a result anticipated, it will be shown that the assumption of bounded acceleration 
effectively eliminates the specific impulse I from consideration in the present analysis. 
This assumption leads to constant acceleration rather than constant-thrust maneuvers 
whose acceleration profiles depend on I. These two are identical in the limit of infinite I, 
and it is shown in reference 4 that I has little effect on AV when maneuvers are com- 
pared on the basis of equal burn time (or equivalently, equal average acceleration). 
Whatever the value of I may be, it is easily shown from equation (2) that the initial and 
average accelerations are related by 


A 


o 


1 



AV 


<A> 


(3) 


This expression may be used to apply the present numerical results to constant-thrust 
finite I maneuvers by identifying the quoted acceleration levels with (A) . Moreover, the 
results apply to both escape and capture maneuvers since A is constant and capture 
maneuvers, after the convention of reference 4, are integrated backward from the desired 
orbital power-off point. 

Equations of motion . - With the use of the foregoing assumptions, a simple system of 
differential equations and boundary conditions may be constructed to describe the trajec- 
tory. In plane polar coordinates, the trajectory variables (illustrated in fig. 2(b)) are the 
polar radius R, central angle 9, velocity V, and path angle a relative to the local hor- 
izontal. The thrust acceleration is described by its magnitude A (assumed bounded) and 

— ► o 

angle of attack (relative to V). The force of gravity has a magnitude p/R and is 

directed toward the origin. These are transformed to the usual dimensionless variables 
by dividing radii, velocities, accelerations, and time, respectively, by the radius, local 
circular velocity, local gravity, and local circular radian period corresponding to the 
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periapse of the intermediate orbit. Polar angles are measured from the periapse ray of 
the major axis. The dimensionless variables are denoted by lower case symbols and 
are related to their physical counterparts by 


r(r) = m 

V 


v(t> = . m. 

V 

c,pa 


a(r) -Mr 

JL 

L R paJ 


V 

T = t _c z pa 

^pa 

G 1 


g = 


JL 

R 2 

ITpaj 


V co = 


Av = 


c, pa 

AV 

7 

c,pa 


(4a) 

(4b) 

(4c) 

(4d) 

(4e) 

(4f) 

(4g) 


These may be used to apply 
the first-order vector state 

the subsequent results to particular cases, 
equations are 

In these units, 


r' = v 

(5a) 


v' = a + g 

(5b) 
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or, in component form, 


r’ = v sin a 


v’ = u 2 a cos Uj - 


sin a 


_ u 2 a sl " U 1 t I'v 2 1 \cos or 
r r 2/ v 


Q, _ V COS O' 


Also, the propulsive effort or Av is defined by 


Av 


■is 


a dr 


(6a) 

(6b) 


(6c) 


(6d) 


(6e) 


This is the criterion to be minimized. The control variables are the thrust angle of 
attack Uj, and the throttle control function Ug. These are selected (as functions of time) 
to yield the lowest possible Av consistent with the boundary conditions to now be de- 
rived. 

Boundary conditions . - The near-planet end of an escape or capture maneuver termi- 
nates in the parking orbit. The boundary conditions appropriate here may be defined in 
terms of the orbit eccentricity ep Q and true anomaly v of the power-on or power-off 
point by the well known conic orbit equations (cf. ref. 11): 
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Since the space vehicle would presumably complete several revolutions in the parking 
orbit, the location of the power-on or power-off point is immaterial from the view- 
point of mission objectives and so can be selected for minimum Av. Thus, equations (6) 
define a set of allowable initial conditions with one degree of freedom - namely, Vp Q . 

In addition, the asymptotic boundary conditions must be met at the sphere of influence. 
That is, for the capture maneuver shown in figure 2(b), the approach hyperbola has a vis- 
viva energy defined by 



where 




A 

v 


pi, 1 


2 

c,pa 


2 


(8a) 


(8b) 


and similarly for escape. Further, a geometric constraint of the form 


0 h + 0 cst “ 2 


(8c) 


must be satisfied by each maneuver. Again, using the conic relations allows the asymp- 
totic polar angle ^cst for the coasting arc to be written as 


9 


cst ^h 


(8d) 


where 


tan m - -hv 

t cO a 


(8e) 


tan 


^h 


hr h r h 

.2 

h - r h 


(8f) 


and (p ^ and cp ^ are, respectively, the asymptotic and power-on or power-off true 
anomalies and h the .’constant) angular momentum of the coasting arc. Equations (8) de- 
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fine a two-degree-of -freedom "target set" of allowable terminal conditions. Its free 
parameters (v and or) were chosen for convenience in applying the transversality condi- 
tions . 


Variational Necessary Conditions 


The optimal controls Uj(t) and UgO") are defined by Pontryagin’s maximum principle, 
which is developed in reference 10. According to this principle, the controls can be op- 
timal only if the Hamiltonian function 


+ ^ 2 V ' + + + ^5 Av ' 


0 ) 


attains its maximum with respect to U^(t) and UgO") for 0 < t < r^. Furthermore, & 
is a constant and its value is zero since in the present case the burning time is not con- 
strained. The state variable derivatives are given by equations (6), and the adjoint 
variables are defined by 


*i- 


r Vr^v r / r 


(10a) 


l//J> = " Wf s * n a + ^9. 


-Ug a sin u^ / ^ 


cos a 


\ r r 2 V 2 / 




cos a 


(10b) 


^ = - 


^v cos a - ^ - -Usin a - ^ 

r 2 \ r r 2 V/ r 


(10c) 


i//^ = - ^ = 0, or i//^ = constant (10d) 

9 0 


^5 = - — = °> or = - 1 
3 Av & 


(10e) 


(The value of -1 is selected as a scale factor since equations (10) are linear and homo- 
geneous in the !/a . The i/a may then be interpreted as the partial derivatives of Av with 
respect to the terminal values of the corresponding state variables.) 
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Optimal control laws. - The Hamiltonian may also be written as 


-# = 


, TT . . sin a , / 2 l\cos a , , V cos a 


rv 


/ sin U 1 \ 

+ (Uga) fi/zg cos Uj + J = F + U2aK (11) 


when equations (6) are substituted into equation (9). Since Uj^t) is not constrained, the 
maximum condition yields 


d# _ o 

au x 



<0 


( 12 ) 


which leads immediately to the optimal steering control law 


tan Uj(t) = 


^ 3 (t) 

^ 2 ( t )v( t ) 


(13) 


On the other hand, the parameter a is interpreted as the maximum acceleration capa- 
bility of the vehicle. Thus, U 2 (r) must be constrained to the range 

0 < U 2 < 1 (14) 


Then since U 2 enters 4€ linearly, its value must switch dis continuously between zero 
and wide open according to the sign of the "thrust switching function" K. That is, 

U 2 = 1 K > 0 ] 

\ (15) 

= 0 K < 0 J 

and is indeterminate in the "singular" case (discussed in appendix A) where K = 0 over 
some finite internal. Thus, aside from the singular case, the optimal trajectories must 
consist of constant-acceleration and coasting arcs. This conclusion, together with the 
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previously mentioned fact that capture maneuvers are integrated backward from the de- 
sired power-off point, means that escape and capture maneuvers are mathematically 
identical; the results which will be shown subsequently apply to both cases. 

As a point of interest it may be noted that by using equations (13) and (lOe), K may be 
expressed as 


K = P - 1 (16a) 

where 



Thus, the vector P with magnitude P and orientation given by equation (13) is seen to 
to the "primer vector” of reference 1, expressed in the present coordinate system. Here 
as in reference 1 the thrust is alined with the primer and is turned on or off according to 
whether the magnitude P is greater or less than unity. 

Trans versality conditions . - Eleven conditions are needed to determine the burning 
time and to specify the simultaneous solutions of equations (6) and (10). Four of these 
are given in terms of v pQ by equations (7), two are given implicitly in equations (8), one 
is given by the choice of scale factor in equation (lOe), and one more is determinable 
from the fact that = 0. The remaining three conditions are obtained by applying the 
transversality conditions at both the initial and final times. (The partial derivatives 
needed in applying this condition are presented in appendix B.) At the orbit end this re- 
sults in 



+ *2 77 + *3 


0 oi 
00 



= 0 


and at the hyperbolic terminal, 





= 0 


(17) 


(18a) 


=0 


(18b) 
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For the case in which 6 ^ is not constrained, these conditions reduce to those presented 
in reference 4. The adjoint variable ^ (corresponding to the polar angle 9) must then 
vanish; that is. 




3Av 



0 


This further implies that ^ = 0 ^ (using eqs. (10c) and (18b) and the relation 

(3r/3v) = -r^v). Thus, in the absence of a geometric constraint, the final thrust direc- 
tion of an optimal escape or capture maneuver must be tangential to at least second order 
in (t^ - t) . This explains in part why the tangential steering law (U^ = 0) has been found 
so close to optimal (in the sense of producing near -minimum values of Av) for the 
optimal-angle maneuvers discussed in reference 4. 


Optimal Trajectory Solutions 

Numerical integration is required to produce solutions in the general finite-thrust 
case, but analytical solutions may also be derived for the limiting cases of impulsive and 
ultra-low thrust. These cases are of interest because of the following: 

(1) They furnish upper and lower bounds for finite-thrust performance. 

(2) The impulsive solution in particular provides good physical insight for a surpris- 
ingly large range of finite-thrust behavior. 

(3) The impulsive results can further be used to generate a simple but effective feed- 
back solution which may be useful for approximate optimization studies. 

These solutions are now presented together with a formal "mean value" result that 
proves useful in interpreting the numerical finite-thrust data. 

Impulsive jthrust. - Reference to figure 2(b) shows that for impulsive thrust the pow- 
ered arc shrinks to a point. In this case, the geometric boundary condition requires that 

0 =O - V = 0 cst = ^=o-^h (19 > 

Thus, 

cos^ - v pQ ) s D = cos (cp^ - <p h ) = cos <p^ cos cp h + sin ^ sin cp h (20) 

Expanding this expression with the help of the conic relations yields 
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( 21 ) 


squaring, and solving the resulting quadratic for p h gives the elements of the coasting 
hyperbola as 



> ( 22 ) 



where 


A = r h (l - D) 

and (23) 

B = r 2 h v 2 Jl - D 2 ) 


The velocity and path angle at the impulse point for both the parking orbit and the coasting 
hyperbola may be computed using equations (7b) and (7c); by the law of cosines the neces- 
sary Av is then 


Av = 




" 2v h v o cos 



(24) 
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and the law of sines yields 


where 


sin U., _ sin U, . _. o 

1, o 1, h sin p 

Vu v„ Av 


(25) 


This is the complete and unique single -impulse solution to the given problem if the initial 
position v is specified in advance. That this solution also satisfies the maximum prin- 
ciple and the terminal transversality condition is now demonstrated. Substituting equa- 
tions (13) into equation (18b) and the result into equation (18a) yields 



1 - 



v^ tan U 


l,h 



^3 



(26) 


(27) 



~^2,h V h tan U l, h 



(28) 


where the partials are defined in appendix B and the scale factor is temporarily 
% 

chosen as = 1. Thus, the terminal transversality conditions are satisfied. Next, the 
previous normalization = -1) may be recovered by computing 



(29) 


5J5 

(obtained by setting K = 0 in eq. (11) for an optimal impulse) and then dividing i// 1 (=1) 
through by -i// g . 

The adjoint variables and j// g are all constant across an impulse, since 

their derivatives are bounded (i.e. , do not depend on a) and the burning time is infini- 
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tesimal. The jump in i/zg across an impulse may be obtained by neglecting all terms not 
dependent on a in equation (10b) and then substituting from equations (6b) and (13): 


K 


IpO o 

2 ~_£(u 2 a sin Uj) = — v' tan Uj 

v 2 v 2 


,2 

^2 v 3 


Hence, 


, . , , 2 dv 

v 3 

The solution of this is 


(30) 


,2 

2 ^3 

+ — = constant (31) 

V 2 

It is clear, when recalling equations (16), that the primer P remains constant across an 
impulse, and if the impulse is optimal with respect to thrust control the constant has unit 
magnitude. In this case i/^ = cos u i> which in equation (31) implies that - v sin Up 
and this is in turn consistent with equation (13). Thus, it is seen that the closed solution 
written previously - that is, equations (19) to (24) - satisfies all the necessary conditions 
of the maximum principle in the case where v is fixed. 

If ^p 0 is left open, it is necessary in addition to satisfy the initial transversality con- 
dition, equation (17). This requires an iterative numerical solution - that is, assuming 
an initial value of v , evaluating equations (19) to (31) and (17), and searching along 
l po until equation (17) is satisfied. In the optimal angle case, 


^4 - ^3 - 0 

Hence, 


^2=1 or Uj = 0 


That is, the impulse is applied tangent to the initial velocity vector. In view of the initial 
transversality condition (eqs. (17) and (B4)), this impulse is to be applied at the periapse. 
The corresponding value of Av is then 
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Av 


min 



The associated optimal value of 6^ is 


opt 


= cos 



(32) 


(33) 


where 


A 


e, = v“ + 1 

n «> 


Ph 


= 1 + e. 


(34) 


Mean value sol ution. - Although numerical methods must be used to obtain the general, 
finite-thrust trajectories, a formal solution for Av can be obtained by converting the 
velocity equation (eq. (6b)) into an energy equation. Multiplying by v and using equa- 
tion (6a) yield 


vv' 



= Ugav cos 


U, 


(35) 


or, recalling equation (8a) , 



= Ugav cos Uj 


Applying the law of the mean to the right side of equation (36) gives 


Av = 


V , - V 

°o, h °o, po 
2(v cos Uj) 


(36) 


(37) 


where 


;v cos Uj) = — f k v cos dT 
T h jT o 


( 38 ) 
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is given by 


and the initial vis -viva energy 




( 39 ) 


Thus, Av may be expressed as the ratio of the energy increment between the initial orbit 
and the asymptotic boundary conditions to the mean path velocity component in the direc- 
tion of thrust. As the energy increment is fixed, it is clear that the control action seeks 
to maximize (v cos Uj) , which is the mean rate of energy addition to the space vehicle. 
The ratio of Av's for finite and impulse thrust, which is termed a "velocity correction 
factor", is thus given by 


f , Av act _< v cosU l>imp 
V Av imp ( v cos u i> act 


(40) 


This factor is often used to account for gravity losses (e. g. , a wide range of data of 
this kind is presented in ref. 4) . Nevertheless, it is evident that control actions and the 
choice of initial conditions will have an explicit effect on f whether a gravitational field 
is present or not. Thus, f in general represents steering losses as much as gravity 
losses. 

Microthrust . - An approximate but useful solution for the low-thrust case is suggested 
by the previous result. It is noted, first, that the energy increment in equation (37) does 
not depend essentially on the shape of the initial orbit. For any elliptic orbit (whose peri- 
apse radius r pa = 1), a circular orbit having the same energy (and period) in absolute 
units has the radius 


r c, equivalent - (^1) 

PO 

The results of reference 4 indicate that the value of (v cos Uj) does not depend strongly 
on e po for the low-thrust spirals. Thus, the familiar "expanding circle" approximation 
now takes the form 

Av microthrust “ v «=,h + O' “ e po) ^ (^) 

(The quantity (1 - e pQ ) 1//2 may be recognized as the circular velocity of the equal-energy 
circular orbit.) This in combination with equation (32) yields 


21 



(43) 


v, microthrust 




(* 


\l/2 
+ e po) 


Note that the result expressed in equation (43) does not depend on the asymptotic direc- 
tion 0 w . By contrast, 0 ^ enters explicitly in the impulsive solution (eqs. (19) to (25)). 
Hence, it may be anticipated that the effect of asymptotic direction constraints will be 
progressively less important as the vehicle acceleration level decreases. 

Variatio nal solutions for finite thrust. - The trajectory and adjoint equations (eqs. (6) 
and (10)), the optimal control laws (eqs. (13) and (15)), and the transversality conditions 
(eqs. (17) and (18)), together with the necessary auxiliary relations, were programmed 
for numerical solution on a digital computer. As previously mentioned, each particular 
case is solved iteratively as a two-point boundary value problem. Each trial integration 
involves 11 subsidiary conditions, of which 6 are stated as physical boundary conditions, 

3 are given in terms of the transversality conditions, 1 is incorporated in the choice of 
scale factor, and the last is determined by the fact that 4€ = 0. Note that some of these 
conditions apply at the initial time and some at the final time. A solution to this problem 
is defined by a complete set of initial conditions which lead, by means of equations (6), 
(10), (13), and (15), to the desired terminal conditions. 

The following procedure is used to accomplish the necessary three-dimensional search 

jJ( 

(using y , Uj q, and as search variables) . Equations (17) and (13) are inverted to 
yield 




2,o 



9l + S -^v n 

30 30 0 


tan U 


1,0 


(44) 


j|( 

when resorting to the normalization = 1. The partials (see appendix B) may be 

computed since v is known. Then 

j-MJ 




o v o 


tan U 


1,0 


(45) 


Finally, if/^ is computed to make # = 0 (see eq. (11)) and is used to renormalize so that 
= -1. With this, all initial values are defined. Integration is stopped when the desired 
value of v^ is attained (see eq. (8a)); this defines the final time t^. A terminal error 
function is constructed to give a measure of error in the remaining conditions; that is, 
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E =£\+l \ + .£g (46) 

where £ ^ and £2 are, respectively, the left sides of equations (18a) and (18b) and (from 
eqs. (8c) and (8d)) 

£ 3 = " [ 0 h + ^h " *iO] 

is the error in asymptotic direction. A multivariable search routine similar to that de- 
scribed in reference 12 is used to drive E to zero. Rapid convergence depends on 
making reasonably good initial guesses for v , Q and 1 fs^. The choice of Q , 

and as search variables results in several advantages. First, they are bounded, 
which means that an exhaustive search could be made. This is important from the con- 
vergence viewpoint and also because no uniqueness theorem is available for the present 
nonlinear type of problem. Secondly, these variables have definite physical meanings, 
which are the same whether finite or impulsive thrust is used. Thus, it is often possible 
to obtain adequate initial guesses by inspecting the easily obtained impulsive results. 

Near-optimal feedback s olutions . - Even when the previous techniques are used, the 
"exact" variational solution technique requires a sophisticated computing facility and sig- 
nificant amounts of machine time. A less exact but faster and more flexible technique is 
developed in appendix C . It relies on a nonoptimal but stable and strongly convergent 
feedback algorithm to satisfy the problem's boundary conditions. Near-minimum Av's 
are then obtained by direct numerical optimization of adjustable parameters in the feed- 
back loop. Comparisons with the variational results indicate that the true minimum pro- 
pulsive velocity increment can be approached within small fractions of 1 percent with 
relative ease. Hence, the feedback technique of appendix C could be used effectively in 
those cases where minor deviations from the true minimum Av can be accepted. 


Matching Conditions for Combined Maneuvers 


The considerations presented above lead to escape or capture maneuvers which yield 
minimum Av for prescribed values of v x and 9^. But to solve the originally stated 
problem, it is also necessary to join pairs of escape and capture maneuvers in such a way 
that (1) the geometric constraint illustrated in figure 2(b) is satisfied, that is, 


V 1 + V 2 = 0 = cos ' 



1 • V ~,2 


v =o, l v °o, 2 


(47) 


23 



and (2) the sum 



+ Av 2 


= minimum 


(48) 


When it is assumed that ? w ^ and v^ 2 have been determined by prior interplanetary 
trajectory calculations, the values of Av depend only on e that is, 



Av. = Av^ j) + Av 2 (0 oO) 2 ) 


(49) 


Thus, it is necessary only that 


d 



SSAv 3SAv rfp _ ^ 

d0 =o, 1 + 77 d0 °°. 2 - 0 


30 


°°, 1 


30 


« 2 


or, using equation (49), 


3Av, 


3Av r 


d0 


~,1 + 


86 «,1 39 «,2 




To satisfy the boundary condition it is necessary that 

d© = J + d0 oo 2 = 0 

J x > 


Hence, equation (50) becomes 


3Av, 


3AVr 


90 oo,i 90 o O) 2 


(50) 


(51) 


(52) 


Now since 0^ is a function of 0^, cp^, and <p^ (see eqs. (8d) to (8f)), Av may be written 
as composite function of those variables; that is, 

Av = A y(9 h ,(p h ,(pj 


and thus 
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9Av _ 9Av + 9Av + 9Av ^ «> 
99 oo 39 h d9 oo d( Ph S0 °° d<p oo "oo 


As is seen from equation (8c), 



= 1 


while from equations (8e) and (8f) 


d( P h d( Po 


'= 0 


90 ^ 3 

OO OO 


Hence, 


9Av 9Av 


90_ 


90. 


= ^4 


and equation (52) reduces to the simple requirement that 

^ 4,1 = ^ 4,2 ( 53 ) 

This condition could also have been obtained by considering the problem as a whole. But 
then it would have been necessary to derive "jump conditions” for the i/a (similar to those 
required for state variable inequality constraints, see ref. 10) to ensure that the desired 
orbit was actually attained. The present approach of matching individually optimal sub- 
trajectories by means of equation (53) not only represents a considerable practical sim- 
plification from the computational viewpoint, but it also produces a range of optimal sin- 
gle trajectories that may be of interest in themselves (e. g. , for a one-way orbiting probe 
mission in which the location of the parking orbit periapse has been prescribed on the 
basis of mission observational requirements) . 


RESULTS AND DISCUSSION 

As previously stated, elliptic parking orbits are currently of interest because they 
offer a large reduction of the propulsive effort or Av when used in place of the more 
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familiar circular orbit. Comparing the values of equation (32) with ep Q = 0 (circular or- 
bit) and ep Q — 1.0 (a highly eccentric ellipse) shows that the saving may approach 41.4 
percent of the reference circular velocity for unconstrained geometry (0 ot open). This 
saving is especially important for trips to the major planets, which have strong gravita- 
tional fields and hence high circular velocities. At Jupiter, for instance, where 
V c « 40 kilometers per second, the AV for a typical capture maneuver (per ref. 7) may 
be reduced from 18 or 19 to 2 or 3 kilometers per second. This will clearly have an 
enormously beneficial effect on vehicle weight. 

On the other hand, the long periods and limited times at low altitude associated with 
highly eccentric orbits may create difficulties in carrying out scientific observations of 
the destination planet. Yet these represent the very purpose and justification for the mis- 
sion. Thus it may be expected that the value of ep Q finally chosen will reflect a com- 
promise between propulsive effort and observational requirements. In figure 3, the Av 



Figure 3. - Relation of parking orbit parameters to mission velocity 
increments. Unconstrained maneuvers; impulsive thrust, 

V R P r L ra- 


sa ving and the orbit apoapse radius are plotted as functions of the orbit period. The 
eccentricity e^ is indicated by markers along the two curves. There is evidently little 
further Av reduction for ep Q >0.9, yet both the period and apoapse radius increase 
very rapidly beyond this point. The value of ep Q = 0.9 therefore appears to be a rea- 
sonable choice and is used for illustrative purposes in the remainder of this discussion. 


26 



Individual Escape or Capture Maneuvers 


The individual maneuvers are of interest in themselves (i.e. , for an orbiting probe 
mission) in addition to being part of the "destination planet problem" for a ^ound trip. 
Unfortunately, Av T s as low as those predicted by equation (32) are not always available 
because of (1) the effect of a constraint on 9 ^ (the optimal value per eq. (33) does not in 
general match the problem geometry) and (2) the AV penalty due to finite thrust (gener- 
ally the propulsive efficiency is reduced when an impulse is replaced by finite thrust). 
These effects are now considered. 

Effect o f asymptotic direction . - Constraints on 0 ^ arise regularly as part qf the 
round-trip problem and may also occur in probe missions if observational requirements 
dictate a particular orbit shape and orientation. The effects of constraining 6 are il- 
lustrated in figure 4 for a relatively low-energy maneuver in which v^ = 0. 25. In fig- 
ure 4(a), the characteristic velocity ratio f (propulsive AV divided by AV m ^ n ) is 
plotted against 9 1 for acceleration levels ranging from impulsive (the lower curve) to 
0.01 local gravity. Clearly, the impulsive limiting solution does represent a lower bound 
of performance for all finite-thrust cases. For impulsive thrust, the minimum value of 
1 for f occurs at 9 ^= 2.8 radians as predicted by equation (33). This case represents 
tangential firing at periapse. A maximum value of f = 5. 8 occurs at 9 ^= 5.2 radians; 
this maneuver consists of firing tangentially at the apoapse. For intermediate cases, the 
impulse does not occur at an apse and is not tangential. Ratioing the values of equa- 
tion (32) for ep Q = 0 and 0.9 shows that an equivalent value of f y for a circular orbit is 
about 7.3. Thus, there is always some saving for the elliptic orbit even when the worst 
possible value of 9 is required. 

A comparison of this impulsive solution (in which f represents steering losses 
only) with the optimal-angle results of reference 4 (where f only includes gravity losses) 
indicates that these two effects can lead to similar performance penalties. In general, 
however, f includes both effects simultaneously, and the particular manner in which 
they combine is a major concern of this report. 

Passing to the finite-thrust cases in figure 4(a), it is seen that dimensionless acceler- 
ation as low as a = 0.01 yield a performance that does not differ markedly from impul- 
sive except in the immediate vicinity of the optimal-angle condition (0 ~ 2.8 rad). 

For a = 0.01, the minimum of 1. 38 for f occurs at the same value of 0^ . Except 
within ±0. 5 radian of this value, the results for a = 0. 01 cannot be distinguished from 
these for impulsive thrust. But when a <0.001, the results (not illustrated here) begin to 
resemble the "microthrust" solution (eqs. (42) and (43)). Thus, for low-thrust systems, 
the results do not so depend strongly on 0 . But for medium- or high-thrust systems 
(a = 0.01), the AV requirement is dominated by 0^ and is relatively insensitive to 
thrust level; that is, gravity losses and steering losses are not directly additive. Hence, 
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90 180 270 360 

Asymptotic central angle, 0^, deg 

(b) Values of adjoint variable ijj 4 = (dAV/dO^). 

Figure 4. - Effects of asymptotic direction constraints on orbit-escape maneuvers. 
Parking orbit eccentricity, ep 0 = 0.9; asymtotic velocity, V^ = 0.25 V c ^ 
constant acceleration. 
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90 180 270 360 

Asymptotic central angle, 0^, deg 


(d) Initial angle of attack. 
Figure 4. - Concluded. 


it appears likely that the effects of constrained 9 ^ could be offset at least partially by 
using a much smaller engine system than would normally be expected. This possibility 
is illustrated later. 

Parts (b), (c), and (d) of figure 4 illustrate the optimal values of the search variables 
discussed previously. In figure 4(b) the adjoint variable is plotted against 9^; it 
simply reflects the slope of the curves shown in figure 4(a) . The initial true anomaly v 
and initial angle of attack Uj g are then shown in figures 4(c) and (d). Because of the 
similarity in shape of these curves, good initial guesses for finite thrust 1 ^ 4 , Uj Q , and 
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v can be obtained from impulsive data, that is, 

^po^<*>)finite thrust ~ v po^ d °° ~ ^impulsive thrust 

where the correlation angle X is about 0. 25 radian for the case illustrated in figure 4(c). 
Values of X appropriate for other values of a and v^ can be obtained by comparing any 
one converged finite-thrust case to the impulsive data. In this fashion a finite-thrust 
"performance map”, such as figure 4, can be generated without major computational dif- 
ficulties. 



L. I I I 

90 180 870 360 

Asymptotic central angle, 0^, deg 

(b) Asymptotic velocity, Voo = 0.5 V C; p a . 

Figure 5. - Comparison of impulsive and finite-thrust maneuvers. 
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Asymptotic central angle, 0^, deg 


(c) Asymptotic velocity, = 1.0 V c pa . 

Figure 5. - Concluded. 

Effect of asymp totic velocity. - Figure 5 illustrates the A V/9^ characteristic for 
values of of 0.25, 0.50, and 1.0, and for accelerations of 0. 001, 0.01, 0.1, and im- 
pulsive. The results are now presented in terms of Av = AV /V (rather than in terms 

c , pa 

of f s (AV/AV . )) so that the Av's for different values of v are made dimension- 

v ' mm" 1:0 

less by the same base number. Aside from this change of base, figure 5(a) represents the 
same low-energy maneuvers (v = 0. 25) that were just described. In comparison with 
this, the higher energy trajectories portrayed in figures 5(b) and (c) display (1) generally 
higher Av levels, (2) less relative variation between the impulsive maxima and minima, 
and (3) a greater separation between the impulsive and low-thrust curves. From this, it 
is concluded that the effects of finite thrust are increasingly more important for high- 
energy maneuvers; that is, the gravity losses are proportionately larger and extend over 
a wider range of 9^. Nevertheless, it is still true that gravity losses and steering losses 
do not add directly together. As previously mentioned, the effects of constraints on 9 ^ 
can be offset to some degree by using a relatively low acceleration and hence a small and 
presumably light engine. 

Characteristics of finite-t hrust trajectories . - In figure 6 are presented some typical 
state and control variable histories for maneuvers where a = 0. 01 and v^ = 0.25. Fig- 
ure 6(a) represents the optimal angle case where 9 is allowed to take on its most ad- 
vantageous value of about 2.8 radians. The state variables r, v, and a, the angle of 
attack control variable Up and the thrust switching function K are plotted against the 
polar angle 9 from v to 9^. The powered arc is indicated by the solid portion of the 
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(a) Unconstrained asymptotic angle (^4 = 0). (b) Constrained asymptotic angle, 0 = 3.34 radians (^4 ~0.2). 

Figure 6 . - Characteristics of finite-thrust trajectories. Parking orbit eccentricity, e po = 0.9; dimensionless acceleration, a = 0.01; 
dimensionless hyperbolic velocity, = 0.25. 


curves while the subsequent coasting arc is dashed. (The dash-dotted curves are dis- 
cussed in appendix C.) 

It is clear that in this case there is very little steering action. That is, is only 
slightly different from zero (tangential thrust) at the beginning and approaches this condi- 
tion ever more closely as the maneuver proceeds. (It was shown in ref. 4 that tangential 
steering is always a close approximation to the optimal when 9 ^ is not constrained.) 
Here cos = 1 throughout the maneuver, but (v) is seen to be appreciably lower than 
for the impulsive case denoted by heavy dots. Hence (recalling eq. (40)), the Av is in- 
creased by some 38 percent. The thrust switching function begins at zero, increases to a 
sizable positive value during the midportion of the trajectory, and then decreases to zero 
again at - that is, just as the prescribed asymptote is attained. The thrust control 
U 2 , which had unit magnitude from t = 0 to r^, is then zero during the final part of the 
trajectory in which the vehicle coasts to the sphere of influence. 

Figure 6(b) presents the same data for a constrained trajectory, that is, 3. 334 
radians = 0.2). Although this is only about 30° beyond the optimal angle, the trajec- 
tory is evidently quite different in all respects . There is a pronounced steering action in 
this case (U^ « -0. 55 rad), which in addition to the late power-on point (about 1.6 rad 
after periapse) provides the extra 30° of turning. This is obtained, however, at the ex- 
pense of doubling the Av as compared to the previous case. That is, with both v and 
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cos Uj low throughout the maneuver, the mean value (v cos Uj) is also low and (recall- 
ing eq. (40)) Av is therefore large. In this case, however, there is little difference in 
Av for the impulsive and finite-thrust cases. 

It is of interest to note the behavior of the thrust switching function K in this case. 
As figure 6(b) shows, its magnitude is very small throughout the entire trajectory. It is 
believed, however, that this is not a singular arc in the mathematical sense discussed in 
appendix A because of the following two reasons : 

(1) Although it cannot be illustrated to the scale of figure 6(b), K in this case has the 
same characteristic arched shape as in figure 6(a) - except scaled down by about three 

_3 

orders of magnitude. The peak value of around 10 is considered significant in compari- 
son to the roundoff errors of numerical integration. The effects of the latter may be 
judged by the fact that which is supposed to be constant and was initially set to zero, 

_ H 

was observed to fluctuate randomly between ±10 . In this case, K contributed about 

IQ' 5 to tt (=F + U 2 aK) and is thus significant by two orders of magnitude. 

(2) In addition to the fact that K behaves qualitatively as it should, there is a ready 

physical explanation for its observed low quantative magnitude: As may be recalled from 
figure 4(a) or 5(a), the Av is quite insensitive to acceleration outside a ±30° band cen- 
tered on That is, for trajectories with strong geometric constraints it simply 

does not matter what thrust is used in the range 0. 01 < a < °°. As may be seen in the 
figure, the powered trajectory takes place at rather high radii (2 < r < 10.3) where the 
gravitational attraction (l/r ) is weak. At a high enough radius, even 0.01 g behaves as an 
impulse, and for an impulse K = 0. (Note that the acceleration ratio in the present ex- 
ample would vary from 0.04 to about 1. 1 if referenced to the local rather than the peri- 
apse gravity field.) 

Thus, it is concluded that this trajectory also is locally optimal with respect to 
thrust control. A study of other cases (not illustrated) suggests that this behavior is quite 
typical of heavily constrained trajectories in general. 


Combined Trajectories 

The main features of individual maneuvers have been indicated, but it remains to 
combine them so as to form an optimal overall trajectory described by the parameters 

V °o, 1’ and V oo, 2 ( reca11 fi S- 2 ( a ))- 

Ma tching procedure . - As shown in the analysis, optimal matching is attained when 
the following two conditions are satisfied: 

9 oo 1 = 0 - d oo 2 ( 55a ) 
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^ 4,1 “ ^ 4,2 


(55b) 


Equations (55) can be satisfied without an iterative search by alining the escape and cap- 
ture 4 against 0 x curves (see fig. 4(b)) in a "back-to-back 71 fashion, so that the sum 
of ordinates is constant and equal to 0. The necessary value of ^4 1 = ^4 2 is then 
read directly from the intersection. This is illustrated in figure 7(a) for an impulsive 
thrust case in which v - = 0.25, v^ 9 = 1.00, and © = 200° or 3. 5 radians (about 80° 



45 60 75 90 105 120 

Asymptotic central angle for escape, 8^ de 9 


(a) Effect of parking orbit orientation. Hyperbolic velocities, 



(b> Comparison of several types of approach and departure 
maneuvers. Parking orbit eccentricity, 0.9; hyperbolic 
velocities, v„ x - v„ 2 = 0. 5 V c _ pa . 

Figure 7. - Combined maneuvers at destination (impulsive AV's). 
Angle between hyperbolic asymptotes, 0 - 3. 5 radians (200°). 
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less than the optimal value defined by eq. (33)). The total propulsive Av^ + Av 2 and the 
adjoint variables y and 2 are Plotted against 0^ ^ and 0 M 2 (note that 0^ ^ 

and 9 m 2 sum to 200° everywhere) over the range from optimal-angle capture to optimal - 
angle escape. The minimum value of Av^ + AVg (point a) and equality of x[/^ ^ and 
«// ^ g dearly occur at the same values of namely 112°. The corresponding 9 X 2 

is then 88°. In general, the optimal distribution of 0 between 0^ ^ and 0^ 2 tends to 
favor the high-energy maneuver; it becomes more nearly even (i.e. , i — 9^ 2) as 
v « 1 * v w 2 ' 1° figure 7(b), the effect of 0 itself is illustrated for pairs of impulsive 
maneuvers where v w ^ = v^ 2 = 0.5. The lower solid curve represents trajectories 
using single -burn maneuvers; every point on this curve was obtained by the previous 
method. The minimum occurs at 0 = 286°, which is compatible with the values of 0^ j 
and 0 ,, («143°) called for by equation (33) in the optimal-angle case. 

o©j ^5 

Comparison of single- and mu ltiple -burn maneuvers . - So far, only single -burn 
maneuvers have been considered. But from references 5 to 7 it is known that alternative 
multiburn maneuvers exist which may also give low Av. Figure 7(b) also presents an 
impulsive Av comparison between optimal single -burn and typical multiburn trajectories. 
In general, the multiburn trajectories consist of optimal-angle escape and capture to- 
gether with one or more auxiliary maneuvers to reorient the ellipse. That is, the mis- 
match angle o = 0 - 0 is accommodated by rotating the ellipse major axis in the or- 
bit plane or by twisting the orbit plane about the major axis. 

For instance, the "circularize-decircularize" maneuver (ref. 5) whose performance 
is indicated by the upper curve, involves two extra maneuvers: (1) an impulse at apoapse 
to attain a circular orbit with coasting to make up the angle ct; and (2) a retro -impulse 
to regain the original eccentricity but with the line of apsides now rotated through the 
angle cr. For this sequence of maneuvers the required auxiliary Av (in addition to the 
optimal-angle escape and capture Av’s) is 

Av aux = 2 ( v c,aa- v aa> (56) 

Lawden's shift maneuver (ref. 1) may be regarded as an optimized version of the 
above, utilizing transfer by an intermediate ellipse rather than a circle. Its performance 
is illustrated by the symmetric broken curve in figure 7(b), the Av shift penalty is given 
approximately by 


Av 


aux 


2 ( v c,aa " v aa> sin 


(57) 


The "apotwist” maneuver (ref. 5) shown by the dotted curve is essentially different in that 
the orbit is not in the plane of the v M 's. It involves only one extra maneuver which occurs 
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at apoapse and serves to rotate the orbit plane about the major axis. In this case, the 
Av penalty is 


av aux * 2v aa sln 


CO 


( 58 ) 


where 


cos u> = 


co S (0 - ») + sin (f>„ , _ opt - |) Sin ^ 2 opt - |) 


COS 


( 0 oo, 1, opt - f) cos ( 0 ~ 2, opt - |) 


This maneuver is applicable only when cr> 0. All of these multi burn maneuvers share the 
feature that Av g ^^ — 0 as e po — 1, for then V c aa and V aa both —0. Thus, the Av 
penalty for nonoptimal 0 can be made arbitrarily small by first using one further im- 
pulse as periapse to attain ep Q « 1. There is then no 0 penalty, and the Av for raising 
ep Q is recovered in the final escape maneuver. Unfortunately, all this requires a great 
deal of time; the global optimum (zero penalty for constrained 0) is attainable only in the 
limit of infinite maneuver time. The data shown in figure 7(b) are for the case where ep Q 
is not changed. Even so, the apse shift maneuvers still involve some time increment be- 
yond that required for single-burn maneuvers. This increment, roughly proportional to cr, 
should be counted as a mission penalty unless useful observations can be made during the 
long, high radius coasts. (The apotwist maneuver has no time penalty if ep Q is not 
changed but can only be used if a > 0.) 

It is noted that in all the preceding cases the Av shift penalty does not depend on v m 
and thus will be comparatively small for very energetic trajectories. Also, for any energy 
level the multiburn maneuvers do seem to offer advantages when 0 is severely con- 
strained. They should be seriously considered whenever a "difficult" value of 0 (e. g. , 
360°) cannot be avoided, especially with high-energy maneuvers. 

On the other hand, the single-burn scheme is always superior in the immediate vi- 
cinity of 0 = 0 j. This is because equation (53) implies 


ao 


©, 


opt 


(59) 


which is deafly not the case for the multiburn maneuvers as they are presently conceived. 
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Sample Application 


The foregoing results can be most readily appreciated by considering a specific mis- 
sion. A good example for this purpose is the 1400-day Jupiter round trip with 200 days 
of stay time illustrated in figure 1 and discussed in reference 7 . This mission is of in- 
terest because it evidently presents the lowest velocity increment sum (SAv) of any 
Jupiter round trip in which a low periapse elliptic parking orbit is used. Although the 
trip time may seem uncomfortably long (a representative Mars or Venus trip takes about 
500 days) , there is a substantial increase in SAv for going to the next lower feasible 
trip time of 1000 days. The 200-day stay at Jupiter would at least provide ample time to 
carry out a significant amount of scientific investigation. 

The planetocentric maneuver data pertaining to this mission were obtained from ref- 
erence 7: 


V 1 = V 0 « 10.00 km/sec 
© = 4. 94 rad 


Since V c ^ « 40 kilometers per second (the reference periapse radius is taken to be 
1. 125 R^), the dimensionless hyperbolic velocities are 


Since the asymptotic speeds are equal and it is assumed that a ca p = a egc , the matching 
condition (eq. (55)) shows immediately that 0 should be divided equally between the 
capture and escape maneuvers. Then 


0 oo i = e oo 2 = 2 ’ 47 rad 

which may be compared to the optimal value of 2.8 radians obtained from equation (33). 
Then, from figure 5 it is seen that AV^ = AVg »0.1 V c ^ or 4.0 kilometers per second, 
at least. Thus, the total required propulsive effort is about 8.0 kilometers per second or 
more depending on the value of a. In comparison, the best available multiburn sequence 
requires 8.6 kilometers per second (using eqs. (32) and (57)) for impulsive thrust, and is 
more severely penalized for gravity losses since it involves two optimal-angle maneuvers. 

The > AV requirement for the capture - orbit - escape sequence is listed in table I 

for several values of a. The values of AV in the first column were obtained 

from the present single-burn maneuvers, while the multiburn AV’s were obtained by 
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TABLE I. - PROPULSIVE VELOCITY INCREMENTS 


AV AT JUPITER 

[Hyperbolic velocities, V^ = V^ g = 1° km/sec; 
angle between hyperbolic asymptotes, 

0 = 4. 94 rad. ] 


Average 

LAV for single 

SAV for multiburn 

vehicle 

burn maneuvers, 

maneuvers using j 

acceleration, 

km/sec 

Lawden's shift, 

(a). 


km/sec 

dimensionless 



oo 

7.92 

8. 57 

0.500 

7.92 

8.57 

.200 

7.92 

8.58 

.100 

7.92 

8.60 

.070 

7.92 


.050 

7.94 


.030 

7.95 


.020 

7.96 

9.23 

.013 

8.07 


.010 

8.35 

10.45 


combining equation (57) with the optimum-angle results presented in reference 4. It may 
be noted that the single-burn maneuvers yield AV's that are lower by 0. 65 to 2. 1 kilo- 
meters per second over the range of a that was considered, with the largest reductions 
at the lowest values of a. 

Representative space-vehicle weights (evaluated on arrival at Jupiter's sphere of in- 
fluence) for these maneuvers are shown in figure 8 as a function of a. Two general types 
of propulsion systems are considered: 

(1) A solid-core nuclear rocket engine, for which a specific impulse of 900 seconds 

and engine thrust to Earth weight ratio of 5:1 are reasonable values 

(2) A hypothetical advanced engine (such as a gas-core nuclear rocket), for which a 

specific impulse of 1800 seconds and a thrust to Earth weight ratio of 1:2 were 
arbitrarily assumed 

For each case it is further assumed that the total payload accelerated away from 
Jupiter is 100 000 kilograms, that no payload is jettisoned at Jupiter, that the vehicle 
structure weight is 25 percent of the initial propellant weight, and that a single propulsive 
stage is used to perform both the capture and escape maneuvers. 

Minimum weights clearly occur at relatively low values of dimensionless (a) , 0. 013 
and 0.025 for the advanced and solid-core engines, respectively. These values of the 
mean acceleration a may be translated into initial accelerations a Q by means of equa- 
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Ratio of average vehicle thrust to weight to Jupiter gravity 
at reference radius of 1. 125 R 4 , <a> 

Figure 8 . - Selection of desirable thrust to weight ratio for escape 
and capture maneuvers at Jupiter. Payload, 100 000 kilograms. 


tion (3). For the solid-core engine this results in an a Q of 0.016 local gravity or 0.033 
Earth surface gravity. In this case, the most appropriate thrust rating is only 13 000 
kilograms or 28 600 pounds. By comparison, an engine four to five times larger would 
be indicated if the multiburn approach is used; this approach would also involve a total 
weight penalty of about 45 000 kilograms or 100 000 pounds at Jupiter arrival or about 
400 000 pounds in the Earth-centered assembly orbit. * 

These results contrast sharply with those that have been obtained previously in the 
case of a low circular parking orbit at Earth, where it has generally been found that the 
most desirable vehicle thrust to weight ratio for a solid-core rocket system is roughly 
0. 3 local gravity, an order of magnitude higher than the present results suggest. Thus, 
the optimal utilization of an elliptic orbit produces two effects tending to decrease the de- 
sired engine size. First there is the reduction in the optimal thrust to weight ratio that 
was just pointed out. There is also the basic weight saving for using an elliptic rather 
than a circular orbit, which (as was illustrated previously in terms of AV) can be quite 
large, especially at the major planets. These effects combined would indicate a reduc- 
tion in the most advantageous engine size by one or possibly two orders of magnitude. 


* Based on a separate Earth escape stage with the previously mentioned performance 
parameters and a Av capability (see ref. 7) of about 8 km/sec. 
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While it would be improper to base sweeping generalizations on this one example, it 
has been plainly demonstrated that the present considerations can have a major impact on 
the question of upper stage engine sizing. A fresh examination of this question, based on 
the methods described herein, seems to be definitely in order. 


CONCLUDING REMARKS 

The problem of transferring between specified hyperbolic asymptotes by way of a 
given planetocentric intermediate orbit has been studied herein. This is a significant 
problem in space trajectory mechanics which presents itself at the destination planet of 
an interplanetary round trip. The very important cases where the intermediate orbit is 
elliptic or hyperbolic have previously been treated only in terms of impulsive thrust; the 
present work is aimed at optimal -control and finite-thrust aspects of this problem. 

The present solutions involve, first, an analysis of minimum effort trajectories for 
escaping from a highly elliptical (or hyperbolic) orbit and attaining a prescribed asymp- 
totic velocity and direction. Analytical solutions to this subproblem are found for the im- 
pulsive and microthrust limits; these provide lower and upper bounds of performance for 
the general, finite -thrust case which is studied numerically. Second, a simple analytical 
"matching condition" is derived which, when satisfied jointly by pairs of escape and cap- 
ture maneuvers, guarantees that the overall capture - orbit - escape sequence is at least 
locally optimal. Numerical results are then presented (in dimensionless form) for repre- 
sentative cases. 

While by no means exhaustive, the present analysis appears to justify several general 
conclusions which have important implications for future space-mission planning and de- 
velopment. 

1. The dramatic AV reductions predicted by prior impulsive studies (refs. 5 to 7) 
for the elliptic parking orbit are maintained and in some cases even enhanced for finite- 
thrust levels as low as 0.01 x (local gravity). Very-low-thrust systems on the other hand 
are penalized by large gravity losses. Thus the elliptic parking orbit is most desirable 
in combination with a medium- to high-thrust propulsion system. 

2. For this thrust regime, the gravity losses due to finite thrust and the "steering 
losses" due to asymptotic direction constraints do not add together directly. That is, the 
phenomenon of gravity loss is confined to optimal- or near-optimal-angle maneuvers . 
Heavily constrained trajectories on the other hand involve significant steering losses but 
there is then little or no additional penalty for finite thrust. 

3. In many of the constrained angle cases, the propulsive effort is very insensitive 
to acceleration. It is shown, however, that the trajectories presented here are actually 
optimal with respect to thrust control. That is, the thrust switching function, though 
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numerically small, is never the less significant compared to roundoff error and exhibits 
the expected type of behavior. The observed insensitivity is easily explained on physical 
grounds when it is noted that the thrusting maneuver invariably occurs at high altitudes 
(where the force of gravity is low) if a strong directional constraint is present. 

4. Single -burn finite -thrust escape and capture maneuvers appear to be superior to 
alternative multiburn finite-thrust maneuvers if the asymptotic direction constraint is not 
excessive. Multiburn solutions yield lower AV for severely constrained trajectories, 
particularly if the asymptotic velocities are high. Further comparisons between the var- 
ious types of maneuvers will depend on accounting for variations in the heliocentric tra- 
jectories, which were assumed fixed for the purposes of this study. It should be noted 
that the problem studied here is actually a subproblem of the overall trajectory. It would 
be very desirable to obtain a matching condition to join adjacent heliocentric and planeto- 
centric arcs in a fashion similar to that used herein for joining escape and capture 
maneuvers. 

5. The preceding conclusions imply that it is feasible and probably desirable to 
employ a small, moderate -thrust engine for the destination -planet maneuvers of future 
space missions. For instance, the example discussed previously indicates that the opti- 
mal thrust to weight ratio for the escape and capture maneuvers of a representative 
round-trip mission to Jupiter (see fig. 1) would be in the range of 0.02 to 0.06 Earth 
gravity. Results such as this could open areas of application for types of propulsion sys- 
tems (such as the gas-core nuclear engine) which might otherwise appear marginal. 
Moreover, the design trade-offs for engines that are even now being developed could be 
affected. The present topic should be examined in much greater detail before commit- 
ments are made for the development of new types of medium- to high-thrust upper stage 
propulsion devices. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 11, 1967, 

789-30-01-01-22. 
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APPENDIX A 


SINGULAR ARCS AND PARAMETER OPTIMIZATION 


As mentioned previously, the thrust control U 2 is indeterminate if the switching 
function K is zero over some finite interval. In this case the loss of information oc- 
casioned by the failure of equation (15) may, in principle, be overcome by noting that the 
condition #€ = 0 now yields two separate relations: 

F = 0 1 

\ (Al) 

K = 0 J 


Since K = 0 over a finite interval, the derivatives K', K", . . . must also vanish over 
this interval. Setting K’ = P* = 0 results in the expression 




(A 2) 


Using equations (6b), (10b), and (10c) allows equation (A2) to be formally solved for the 
throttle control U 2 in the form 


U 2 a 


sin U 1 \J/n 

iZ/ 0 cos U. 

v 2 v 2 ] 


= X(r,v, ^ 4 ) 


(A 3) 


Unfortunately, the bracketed term vanishes identically whenever the optimal steering law 
of equation (13) is used. Although it would be possible to continue in this manner (i. e. , 
by setting X' =0), the resulting expression is so cumbersome as to be completely im- 
practical for either calculation or interpretation. Therefore, an alternate approach is 
proposed here. Since this report is primarily concerned with the class of single-burn 
maneuvers, the quantity U 2 a = b taken together is treated as a design parameter, the 
constant value of which is to be set before thrusting begins and cannot be changed there- 
after. The parameter b is now regarded as a new state variable defined by 

b* = 0 (A4a) 


b(0) = b(T h ) , unconstrained 


(A 4b) 
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Corresponding to b there is a new adjoint variable xp-^ defined by 

xpL = - — = 1 - P = -K (A 5a) 

D ab 

^ b (0) = ^ b (r h ) = 0 (A 5b) 

where the boundary values (eq. (A 5b)) result from application of the transversality condi- 
tion. 

Evidently, equations (A 5) are satisfied for any arc that is singular from beginning to 
end. That is, any constant level of acceleration that produces a singular arc is optimal 
when viewed as a design parameter. 

Equations (A 5) are also satisfied in the case of impulsive thrust, since then the burn- 
ing time is zero. It is also of interest to note that if Ug were retained as an active con- 


trol, the conditions for optimal parameter a are 

a' = 0 (A 6a) 

a(0) =a(Tj 1 ), unconstrained (A6b) 

K = U 2^ - p ) (A7a) 

* a (0) = ^ a (T h ) = 0 (A 7 b) 

In this case, the quantity UgCl - P) is nonpositive (since Ug = 0 if P = 1) and hence, 

^ a (T h ) ^ 0 (A 8) 


provided that ipAO) = 0. Thus, if unlimited throttling control is admitted, the only pos- 

a 

sible form of optimal trajectory in the present context consists of either impulses or 
singular arcs in combination with coasting (zero thrust) arcs. That is, no finite-thrust 
level can be optimal except as it occurs in a singular arc. 

In other cases, finite acceleration produces a negative value of ^(t^) which implies 
that Av could have been reduced by raising a. This is a general conclusion and is not 
limited to the present problem which after all is distinguished from other space trajec- 
tory problems (such as orbit transfer) only by the form of the boundary and transver- 
sality conditions. 
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APPENDIX B 


PARTIAL DERIVATIVES FOR THE TRANSVERSALITY CONDITIONS 


The transversality condition at r Q requires that 


xb • = 0 

r o o 


(Bl) 


where T Q is the tangent vector to the inner manifold which is defined (after substituting 
eq. (7a) into eqs. (7b) to (7e)) by the following relations: 


1 + e, 


r = 


po 


° 1 - c po cos V 


V ■ | + 2e P° C0S V ■ 1 , e 


0 o y po 


1/2 


PO 


, e po sin y po 

tan a = - — 

° 1 + e cos v 

po po 


Av o = 0 > 


(B2) 


Now the components of T Q are 


9r o 9v o 9 % 90 o 


9y po 9v po 9y po 9 > 


, 0 


(B3) 


Carrying out the indicated differentiations gives equation (Bl) in the form 


*10 7T~ (e po sln v + *20 

X + 6 „ 


po 


-1 


Jo (1 + e po>J 


(e nn sin v ) 
' po po' 


^30 


2 2 

e r cos a ^ 

— ~ ( COS + O 


V + W 


po po 7 


+ ^,=0 (B4) 
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which may be written after simplification as 


sin v = 


-1 


, 2 ^20 
^10 r o " 


cos v + e \ 1 + e nn 

^30 r o cos “ “o 1 , P ° P + 


1 + e. 


po 


“PO 


(B5) 


At time t^, the hyperbolic manifold is defined using v and a as surface coordinates: 


r = r(v) 


v = v 


a = a 




> 


9 = 0(v, a) 

In this case the transversality condition requires that 

■\ 

?h • ? lh = 0 
^h ' ^2h = 0 

where the two tangent vectors are 


0r 


00 


rp _ 'U1 I fv UU 

A lh " l— » L > u > — 


0V 0V 




*2h “ (' 


o, o, i, 

da 


Now from equation (8a) , 


(B6) 


(B7) 


(B8) 


r h = 


2 2 
v h " v °° 


hence, 
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9r 

3v 


■ 


(B9) 


Also, from equation (8c) 


0 h = 0 co+<P h " 


where 


and (8f) 

with respect 

to V 

90 h_ 

dcp h _ 

d( P oo 


1 

XJ 

> 

CO 

dv h 

0v h 

> 

30 h _ 

d<P h _ 

d( P OO 


ea h 

9 “h 

9 “h 

J 


(BIO) 


d( P h 


sin 2 <p h 


0v h sin a h cos “h\v 2 


dcp 


da 


- = sin 2 <p h |cot <p h (cot Q? h - tan af h ) + 2 


3<PoO 81112 ^OO / v oo l 


3v h 2 cos * h Voo/ 


d< Pc 

da. 


= -sin cp^ cos (p^t an o- h 


(Bll) 


These are the relations needed for the evaluation of equations (17) and (18). 
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APPENDIX C 


NEAR-OPTIMAL FEEDBACK SOLUTION 

The previously described technique of solution, consisting of a numerical search to 
zero the variational "terminal error function" (see eq. (46)), is quite appropriate for the 
purposes discussed in the text because it precisely defines the minimal Av and other 
features of the optimal trajectories. Unfortunately, it requires extensive computational 
effort, and there are numerous cases in which a less precise but faster and more flexible 
technique of solution would be more suitable. 

One such technique is presented in this appendix. A suboptimal but strongly conver- 
gent feedback algorithm is used to compute the steering signal U^(r). This provides a 
feasible solution of the two-point boundary value problem in every case. That is, the de- 
sired terminal values are always attained (though not for minimum Av) in spite of initial- 
condition errors or other disturbances. Optimality is approached by introducing an arbi- 
trary power series into the feedback loop to temporarily bias the aiming data (9^) fed into 
the steering algorithm. The coefficients of the series are then determined by direct nu- 
merical search to yield minimum Av. Thus, the variational two-point boundary value 
problem is eliminated in favor of a more tractable and conceptually simpler problem in 
the ordinary calculus. As will be seen, the "true" (i.e., variational) minimum Av can 
be attained within practical tolerances by considering only the first few terms of the bias- 
ing series. 


A Basic Feedback Algorithm 

To eliminate the difficulties of the two-point boundary value problem, an alternative 
scheme must necessarily rely on the inherent goal-seeking nature of feedback control. 

The impulsive solution described previously does in fact possess the necessary closed- 
loop or goal-seeking character. That is, the proper angle of attack may be found by 
means of equations (19) to (25) as a function of the current values of the state variables 
(r, v, a, and 9) and of the desired terminal conditions (v^ and 9^). 

This scheme may be extended immediately to the case of finite thrust. At every in- 
stant, the "velocity to gain" vector Av^ whose components are given by equations (24) 
and (25) is computed just as though the remainder of the trajectory were to be performed 
impulsively. Although Av^ must now be computed repeatedly during the maneuver, each 
calculation requires no information beyond the current state variables and the desired 
terminal conditions. This is, in fact, a feedback control policy for the stated problem; 
the computational algorithm is given by equations (19) to (25) with r^, v Q , <y , and v 


47 



replaced by the current values of the state variables r(r), v(t), or(r), and 9(t). 

This is sometimes referred to as the "cross-product” rule, because (in 3 dimen- 
sions) the alinement of a and Av corresponds to the condition 

ax Av = 0 (Cl) 

O 

As was shown in the ANALYSIS, this policy is actually optimal in the impulsive limit; it 
is also known to be fairly efficient for high-thrust systems such as chemical rockets (see 
ref. 13). Moreover, it produces the largest possible instantaneous rate of decrease of 
AVg, regardless of the acceleration level. This is easily seen by taking the time deriva- 
tive of equation (24): 

AVg(AVg)’ = v^v^ - v cos 0) + v'(v - v^ cos 0) + sin 0(0;^ - <*') (C2) 

In equation (6) , it is noted that only v* and o'* depend explicitly on Uj. Thus, 

i = 0 = (v - v. cos 0) - (vv. sin 0) 

aUj aUj n 3U X n 

or 


w. sin 0 -a sin U., 

= = -v tan Uj 

v - v^ cos 0 a cos Uj 

v 


3v* 



da' 


3Uj 


Hence, 


tan U. 


v^ sin 0 

max v^ cos 0 - v 

Av’ 

g 


That this represents a maximum rate of decrease is seen by noting that 


3 2 Av’ 

s = a cos U., 

9 1 

3UT 


(v^ cos 0 - v) + 


v 2 sin 2 0 


V h cos 0 - V 


(C3) 


(C4) 
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Since s v, cos and the quantity v^ cos j3 - v has the same sign. 

Convergence and Stability 

For the same reason, the convergence of this algorithm is assured for sufficiently 
high thrust. In that case, equation (C2) is dominated by the terms 

( a sin U,\ 

-j (C5) 

and AVg approaches zero monotonically. Even when the terms in equation (C5) are not 
clearly dominant, it appears that convergence can be attained without difficulty as long as 
the steering signal remains within the first or fourth quadrants. That is, it was shown in 
the Mea n value solution section (p. 20) that the instantaneous rate of adding energy per 
unit mass is 



= av cos Uj 


(C6) 


2 

and hence v m (t) increases monotonically toward its desired value so long as cos is 
positive (U^ in quadrant 1 or 4). After local escape energy (v^ = 0) has been attained in 
this manner, the vehicle must necessarily begin to traverse an arc along which r in- 
creases monotonically, which implies that r’ > 0. From equation (8a) it is recalled that 



2 

°°, des 



(C7) 


so that 


I 


V 


h " 



(C8) 


Hence, when r' is positive and Uj is in quadrant 1 or 4, the velocity to gain Av^ 
approaches zero monotonically regardless of thrust level, because every term in equa- 
tion (C2) is then negative. (Actually, it is clear that these conditions need only be satis- 
fied over the terminal portion of the trajectory to guarantee convergence.) As a point of 
interest, it is clear from the preceding discussion that (Av ) may serve as a Lyapounov 

o 
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stability function for this problem (see ref. 14) . That is, a region of asymptotic stability 
exists for sufficiently large r, and this region is attainable by every steering control 
which lies entirely in the first and/or fourth quadrants. 


Optimal Compensation 

Although the preceding "cross-product” algorithm is optimal only in the limit of very 
high (impulsive) thrust, its strong convergence properties suggest that it would also be 
useful for the medium -thrust regime of interest here. Unfortunately, initial numerical 
simulations proved disappointing. Although the stated convergence properties were 
clearly demonstrated, a serious reduction in propulsive efficiency was noted for acceler- 
ations lower than about 0. 2 local gravity. The reason for this is not hard to find and can 
easily be understood by considering figure 9. There, three impulsive trajectories are 
illustrated which depart from different points on an ellipse and attain the optimal value of 
9 m defined by equation (33). As has been previously noted, the minimum energy trajec- 
tory departs from periapse (position 1) and employs tangential steering, = 0. If the 
impulse occurs before periapse (position 2) a sizable negative angle of attack is required 
to produce a trajectory that will attain the necessary value of 0 . Similarly, a late im- 
pulse (point 3) requires a positive angle of attack. Now it was shown in reference 4 that 
a finite-thrust arc is generally "centered" on the location of the equivalent impulse. This 
means that for an optimal-angle trajectory, the powered arc is distributed into roughly 



Figure 9. - Impulsive trajectories attaining the same asymptotic direction. 
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equal central angles ahead of and behind the original periapse. For vehicle accelerations 
of the order of 0.01 local gravity the "lead angle" v is of the order of -1 to -2 radians. 
Thus, initially the elementary controller would command a large negative angle of attack 
(« -1.4 rad); the commanded value would gradually approach zero as the original periapse 
position is approached and finally increase to a sizable positive value. In view of equa- 
tion (38), this behavior, illustrated in figure 10(a) for a = 0.02, leads to a large increase 
in Av. 

On the other hand, the preceding algorithm (eqs. (19) to (25)) can be inverted to com- 
pute instantaneous values of 9 m along an arbitrary trajectory. When this procedure is 
applied to the optimal -angle trajectories presented in reference 4, it is found that, in the 
thrust regime of interest, the optimal value of 9 ^ is initially biased by -0. 1 to -0.4 
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(c) Optimum bias signal. 


Figure 10. - Comparison of controi histories. Dimensionless acceleration, 
a = 0.02; parking orbit eccentricity, e p g * 0.9; asymptotic central angle, 
8^ = 2.8 radians; hyperbolic excess velocity, v^ = 0.25. 
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radians and increases asymptotically toward the desired value. A typical example is il- 
lustrated in figure 10(b), where ^ is plotted against t for an optimal trajectory 
with e po = 0.9, a = 0.02, v w = 0.25, and 9 M dgg = 2.8 radians. The bias signal 


b d °°, des 0 «>, opt 


(C9) 


is plotted against Av g in figure 10(c). It is clear from this that the optimal bias signal 
can be represented in the form 


b opt (T) = £ V" 

n=l 


(CIO) 


where 


Av (t) 

Y = — i— 

AvjO) 

b 

and can in fact be rather well approximated by keeping only the first few terms of this 
series. Equation (CIO) is used as the basis of a compensator. This modifies the 9 x 
signal fed into the cross product algorithm, equations (19) to (25), in such a way as to 
produce a value of Uj(t) which, in the absence of perturbations, is very nearly the same 
as U 1 (r) 0 p t . That is, 


des + b opt< T > < CU > 

where b Qpt is defined by equation (CIO). It is postulated that optimal values of the b^ 
determined in a nominal unperturbed case by the Ritz method of numerical optimization 
(see ref. 15), will remain very nearly optimal when a small trajectory perturbation is 
introduced. Example perturbations include (1) errors in initial conditions, (2) off-nominal 
engine performance, and (3) errors in the form of the state equations (eqs. (6)). The 
present, compensated cross-product algorithm is shown to be capable of withstanding 
rather large perturbations of these kinds without violating the desired boundary values 
(v TO deg and 9 ^ deg ) and with only a small penalty as compared to the variational solu- 
tion in which the same perturbation is accounted for. 
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Comparison of Open-Loop and Closed-Loop Results 


The first question to be discussed is whether, or how fast, the Ritz compensated 
cross-product algorithm can be made to approach the corresponding open-loop optimal 
results presented in the text. Toward this end, consider again the optimal-angle escape 
problem for e p O = 0*9 and v oo =0.25. (By eq. (33), the value of 9 ^ Q pj. = 2. 8 rad.) In 
figure 11, the Av penalty for using the present controller rather than open-loop varia- 



Figure 11. - Effect of truncating compensator 
series. Parking orbit eccentricity, e pq = 0.9; 
asymptotic central angle, 6^, ■ 2. 8 radians; 
hyperbolic excess velocity, = 0.25. 


tional steering is presented as a function of N, the index of the compensator series. 
Values of a of 0. 1, 0.02, and 0.01 local gravity are represented by the lower, middle, 
and upper curves, respectively. For relatively high acceleration, the Av penalty is 
small even at N = 0, and for practical purposes it is eliminated by including one term in 
the series. For a = 0.02, the initial penalty is larger and two terms are needed to sen- 
sibly eliminate it. (Its control history was compared with that for N = 0 and variational 
steering in fig. 10.) The case where a = 0.01 is considerably more difficult, requiring 
at least N = 3 to bring the penalty down to an acceptable level. 

From these observations, the following are immediately inferred: 

(1) For a given value of N, the penalty decreases rapidly as a increases. (Recall 
that the uncompensated policy was shown to be optimal for impulsive thrust.) 

(2) For a given value of a, the Av penalty is reduced by nearly an order of magni- 
tude when N is increased by one unit. 

(3) In the absence of perturbations, the variational Av can be approached within ar- 
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bitrary tolerances by choosing a large enough value of N", if the coefficients b^ . . . bjj 
are selected optimally. 

Some very desirable properties of the compensated policy were demonstrated pre- 
viously for the case where 9 m is not constrained. These advantages also remain in full 
force when angular constraints are imposed (i.e. , values of 9^ different from that given 
by eq. (33)). Refer again to figure 5 where values of the propulsive Av, normalized by 
equation (32), are plotted as functions of 9 ^ for the entire range 0 < 0^ :£ 2tt. The 
lower, solid curves represent the optimal impulsive solution discussed previously; the 
others represent various finite values of a. For fixed values of a as low as 0.01, the 
present closed-loop controller yields Av's which are virtually as low as those obtained 
by variational means. That is, with the compensator series N < 3, the Av penalties 
are not discernible to the scale of figure 5. 

It may also be recalled that, for a ^ 0.01, the low-thrust trajectories are virtually 
as efficient as the impulsive ones except in the immediate vicinity of 9 ^ In the 
nonoptimal- angle region, the results are quite insensitive to the acceleration level, so 
much so that 0.01 local ”g" behaves almost like an impulse. In this insensitive region, 
the Av’s resulting from the present controller approach the variational minimum even 
more rapidly than in the optimal-angle case. For example, when a = 0.01 and 
0^ = 3. 34 radians, the variational Av was approached within 0. 3 percent by using only 
one term in the compensator series. (For the optimal -angle case, three terms had to be 
used in achieve a comparable result.) In this sense the present scheme appears even 
more advantageous in the nonoptimal -angle region than it did for optimal -angle trajec- 
tories. 

The two trajectories just cited are the ones that were illustrated geometrically in 
figure 6. Figure 6(a) represents the optimal-angle case with 0^ = 2.8, and figure 6(b) 
covers the trajectory where 9 00 = 3. 34. Recall that r,v, a , and the control variables are 
plotted as functions of the central angle 9. Solid and dashed lines represent the varia- 
tional open-loop results, and dotted lines denote the closed-loop results with N = 2. In 
either case, the only difference discernible to the scale of figure 6 is in the steering con- 
trol history U^. As has been previously noted, the present control action oscillates 
slowly about the variational control history. In the optimal-angle case, the oscillations 
occur around a near-zero nominal value and hence produce a definite second-order effect 
on (v cos U^) . For the constrained case, the oscillations take place around a strong 
nominal steering action; the effect is hence to increase cos Uj in one place and de- 
crease it in another. As a net result there is a cancellation effect, and the change in 
(v cos Uj) is considerably smaller than second order. 
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Effect of Perturbations 


The previous results confirm the analytical predictions that were made about the ef- 
ficiency of the Ritz -compensated cross-product algorithm. It finally remains to demon- 
strate its closed-loop aspect - that is, its ability to attain the prescribed boundary con- 
ditions in spite of perturbations (such as a poor choice of initial search variables) which 
would cause the open-loop variational procedure to diverge. Toward this end, the nomi- 
nal trajectory depicted in figure 6(a) (optimal-angle case) was subjected to sizable varia- 
tions in the initial conditions and in the acceleration level a. The results are compared 
with corresponding variational runs in table n. For both schemes, the propulsive Av 


TABLE n. - EFFECT OF TRAJECTORY PERTURBATIONS ON GRAVITY AND 
STEERING LOSS PROPULSIVE VELOCITY INCREMENT 
CORRECTION FACTOR f y 

[Parking orbit eccentricity, e pQ , 0.9; dimensionless hyperbolic velocity, v oo = 0.25; 
asymptotic direction, 0^ = 2.80 rad; dimensionless acceleration, a = 0.02.] 


Perturbation 

Present results (for 

Variational results 

Difference, 


N = 2; constant 

(reiterated) 

percent 


b^ and b 2 ) 



None 

1.1391 

1.1373 

0.16 

Initial true anomaly, 

1.1764 

1.1690 

.63 

c o 

v : 6 

po 

Initial true anomaly, 

1.1589 

1.1504 

.74 

e o 

v : -6 

po 

Dimensionless acceleration, 

1.1911 

1.1842 

.58 

a: -10 percent 





(normalized by eq. (32)) is presented for the nominal case and with perturbations of ±6° 
on the initial position Pp Q and -10 percent on the acceleration level. In all cases, the 
present algorithm, using the same values of bj and bg determined for the nominal 
case, had no difficulty in attaining the required boundary values, and involved only frac- 
tional percentage penalties in Av as compared to the reiterated variational trajectories. 
As was pointed out previously, reiteration is required in the variational case to recover 
the desired boundary values; this implies a priori knowledge of the perturbation. It was 
noted that the Av difference between the present and variational results could be reduced 
to 0.16 percent or less in the perturbed cases by reoptimizing bj and bg; this also as- 
sumes a priori knowledge of the perturbation and is perhaps a more representative com- 
parison than that indicated in table n . 
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